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ABSTRACT 

Context. The soft X-ray emission from high density plasma observed in several CTTS is usually associated with the accretion process. 
However, it is still unclear whether this high density "cool" plasma is heated in the accretion shock, or if it is coronal plasma fed or 
modified by the accretion process. 

Aims. We conducted a coordinated quasi-simultaneous optical and X-ray observing campaign of the CTTS V2129 Oph. In this 
paper we analyze Chandra grating spectrometer data and aim at correlating the observed X-ray emitting plasma components with the 
characteristics of the accretion process and of the stellar magnetic field constrained through simultaneous optical observations. 
Methods. We analyze a 200 ks Chandra/HETGS observation, subdivided into two 100 ks segments, of the CTTS V2129 Oph. For the 
two observing segments, corresponding to two different phases within one stellar rotation, we measure the density of the cool plasma 
component and the emission measure distribution. 

Results. The X-ray emitting plasma covers a wide range of temperatures: from 2 up to 34 MK. The cool plasma component of 
V2129 Oph (r a; 3 - 4MK) varies between the two segments of the Chandra observation: high density plasma (logA'c = 12.1^j j) 
with high EM at ~ 3-4 MK is present during the first observing segment; during the second segment this plasma component has 
lower EM and lower density (log We < 11.5), although the statistical significance of these differences is marginal. Hotter plasma 
components, T > 10 MK, show variability on short time scales (~ 10 ks), typical of coronal plasma. A clear flare, detected during the 
first segment, could be located in a large coronal loop (half length > 3 ^,,). 

Conclusions. Our observation provides further confirmation that the dense cool plasma at a few MK in CTTS is material heated in 
the accretion shock. The variability of this cool plasma component on V2129 Oph may be explained in terms of X-rays emitted in the 
accretion shock and seen with different viewing angles at the two rotational phases probed by our observation. In particular, during 
the first time interval direct view of the shock region is possible, while, during the second, the accretion funnel itself intersects the 
line of sight to the shock region, preventing us from observing the accretion-driven X-rays. 

Key words. Stars: circumstellar matter - Stars: coronae - Stars: individual: V2129 Oph - Stars: pre-main sequence - Stars: variables: 
T Tauri - X-rays: stars 



1. Introduction 



During the pre-main sequence phase low mass stars evolve from 
the embedded protostellar stage to the main sequence. In the first 
few million year of this transformation they undergo a key phase 
in which, while gravitationally contracting, they accrete mate- 
rial from a circumstellar disk (see e.g. |Bertout||1989[ Kenyon 
& Hartmannj 1995 1. During this phase low mass stars are named 
Classical T Tauri stars (CTTS). When the accretion process ends, 
stars, still in the pre-main sequence but lacking the strong line 
emission characteristic of the previous phase, are classified as 
weak line T Tauri stars (WTTS). 



Recently it has been evidenced that the accretion process in 
CTTS is associated with the presence of high density plasma at 
temperature of a few MK radiating soft X-rays. Many fundamen- 
tal aspects of this link are still debated. Moreover understanding 
the accretion process and its high-energy radiation is a key step 
in the understanding of the physics of pre-main sequence low 
mass stars and of their circumstellar environment, because ac- 
cretion strongly affects several aspects of the central star itself: 
accretion determines stellar rotationa l velocities (jBouvier et al.| 
11993), and affects coronal activity ( jPreibisch et al. 20051); the 
accretion history in the early stages of stellar evolution might 
modify the stellar radius, and hence its luminosity, up to ages of 
a few Myr ( ,Bai'affe et al.,2009| ); in addition the accretion process, 
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because of the intense ionizing/heating radiation produced (UV 
and X-rays), influences the physics and the chemistry of circum- 
stellar disks where planets are forming and may even determine 
their lifetime (e.g. [Alexander et a l. 2004, Ercolano et al.|20(j9] l. 

The accretion process is a complex phenomenon in which a 
fundamental role is played by the stellar magnetic field: it trun- 
cates the inner part of the circumstellar disk and guides infalling 
material along its flux tubes toward the central star ( Konigl 
1991[ l. The magnetic field, interacting with the inner disk and 
channeling both inflowing (i.e. accreting) and outflowing mate- 
rial, mediates the angular momentum exchange between the star 
and the circumstellar material (S hu et al.|1 994). 

The accreting material from the inner disk approaches the 
star in (almost) free fall. The impact with the dense stellar at- 
mosphere produces shocks at the foot-points of the accretion 
streams. Most of the kinetic energy of the infalling material is 
thus converted into heat, generating a region at higher tempera- 
ture with respect to the surrounding photosphere, a so-called hot 
spot. Evidence of such hot spots atT ~ 10 000 K is provided by 
optical/UV photometry and spectroscopy of CTTS. 

From a theoretical point of view accretion onto CTTS can 
produce significant amounts of X-rays, since the accreting ma- 
terial should impact onto the stellar surface with a high infall 
velocity (v ~ 300 - 500 km s"'), being thus heated in the accre- 
tion shocks up to temperatures of ~ 2 - 6 MK. Such hot mate- 
rial would radiate mostly in the soft X-ray band (0.1-1 .0 ke V, 
GuUbring|1994 1. The predicted X-ray luminosities in this band 
are large even for moderate mass ac cretion rates (i.e. Lx 



Sacco et al. 



20081. The 



10^° ergs-' for M ~ lO^'^Moyf 
observed Lx from CTTS in the 0.1-1.0 keV band are, however, 
significantly lower than these predictions. 

Observationally CTTS are X-ray bright when one considers 
a broader 0.1 - lOkeV band due to their very active coronae. 
In fact their intense magnetic fields, other than controlling ac- 
cretion, are also responsible for the heating and confinement of 
large amounts of coronal plasma. Most of the X-ray emission 
from CTTS is actually due to very hot plasmas (T > 10 MK, 
much higher than that produced in accretion shocks) and is char- 
acterized by frequent flaring activity, clear evidence of coronal 
emission. Measurements of the density of the "cool" plasma 
(T ~ 3 - 4MK) in CTTS are, however, at odds with a coro- 
nal origin of at least part of the X-ray emission: in 6 out of the 7 
CTTS for which plasma densities have been measured ( Kastner' 
et al. |2002tlStelzer & Schmitt 2004; Schmitt et al. 2005 , Giinther 
et al . 2006 ; Argiroffi et al. 2007; Gudel et al. 2007; Robrade & 
Schmitt(2007;|Huenemoerder et al.)2007| ) the cool plasma is sig- 
nificantly denser (A^e > 10^'cm~^), than ever observed in the 



same-temperature plasma of older late type stars ( jTesta et al. 



|2004j) or in disk-be aring young stars without evidence of accre- 
tion ( [Kastner et al. 2004 ; Gi inther et al.|2010) . Such high density 
plasma cannot be interpreted in terms of steady coronal loops 
( [Argiroffi et al.|2009] ). These findings indicate that the high den- 
sity plasma component in CTTS is associated with the accretion 
process. In particular it is suggested that this plasma is material 
heated in the accretion shock. 

In addition to these studies focused on the plasma densities, 
interesting results have also been obtained by Giidel & Telleschi[ 



( [2007 1 that, investigating soft (~ 0.5 keV) X-ray flux level of 
CTTS and WTTS, proved that CTTS have higher soft X-ray 
fluxes than WTTS. They suggest that this soft excess is due to 
shock heated plasma, or to magnetically confined plasma inter- 
acting with (and/or originating from) accretion funnels. 



optical luminosities are correlated in CTTS (and not correlated 
in WTTS). Such a correlation may be explained if X-rays and 
optical radiation are both produced in the accretion shock and are 
modulated by a time-variable accretion rate and/or by the stel- 



lar rotation. Supporting an idea previously proposed by Gregory 
| et al.[ ( |2007] ), [Flaccomio et aL] however, favored a different sce- 
nario in which both the optical and soft X-ray variability are 
produced by a variable absorption, likely due to disk warps, as 
proposed by [Bouvier e"taL| ( |2D07) for the CTTS AA Tau. 

All these results clearly indicate that the soft X-ray emis- 
sion from high density plasma and the accretion process are 
somehow linked. Different scenarios have been depicted to ex- 
plain this link, but many fundamental points are still discussed. 
The final evidence, clarifying whether the high density cool 
plasma component is material heated in the accretion shock 
(Kastner et al.^2002) , or whether it is coronal plasma somehow 



modified/affected by the accretion process (Giidel et al. 2007 
Brickhouse et al. 2010| l, is still lacking. Moreover a fundamental 



role is possibly played by a variable or non-uniform absorption 
produced by accretion streams and inner disk warps in CTTS. 

Very recently a few high quality X-ray datasets of CTTS pro- 
vided tighter observational constraints and evidenced important 
discrepancies with the predictions of current accretion-shock 
models. The high density cool plasma components of TW Hya 
cannot be explained by the shock of a single accretion stream 
( Brickhouse et al.|[2010[), b ut models with multiple streams are 
required ( Sacco et al.|2010 l. Moreover the shock-heated plasma 
model, applied to the observed X-rays from CTTS, provides sig- 
nificantly underestimated accretion rates with respect to mea- 
surements based on optical indicators (a factor of 10 or even 
more, Argiroffi et al. 2009 Curran et al. 201 1 1. Local absorption 
of X-rays by the surrounding stellar atmosphere (Drake 2005 ) or 
by the pre-shock infalling material could play a fundamental role 
in determining the amount of X-ray emission escaping from the 
shock region ( Sacco et al.|2010| l. 

The determination of the X-ray luminosity that emerges from 
the shock region, other than for interpreting the observations, is 
fundamental to: /) constrain the amount of high-energy radiation 
that heats and ionizes the circumstellar environment, and, //) de- 
lineate the energy balance in the shock region. 

In order to address the above issues, and hence contribute to 
the understanding of the accretion phenomenon as a whole, we 
have planned and performed a campaign of simultaneous mul- 
tiwavelength observation of a CTTS. The immediate objective 
is to simultaneously constrain and correlate the characteristics 
of the accretion process, of the stellar magnetic field, and of 
the X-ray emitting plasma, including both the coronal and the 
accretion-related components (Gregor y et al.|20()9] l. 

We selected for this study V2129 Oph, a CTTS in the 
p Ophiuchi star forming cloud and already one of the targets 
of the Magnetic Protostars and Planets (MaPP) project ( |Donati[ 
et al.[[20() 9'), whose main goal is studying the large-scale mag- 



Flaccomio et al. ( 2010| l, by observing the NGC 2264 star 



forming region, found that time variability of the soft X-ray and 



netic fields of CTTS using time-series of polarized Zeeman sig- 
natures. 

In our study of V2129 Oph we benefit from several datasets, 
acquired simultaneously or quasi-simultaneously. V2129 Oph 
was monitored by optical spectroscopy and photometry from 
early May 2009 to mid July 2009. The observation set includes: 
high resolution X-ray spectroscopy obtained with the Chandra 
High Energy Transmission Grating Spectrometer (HETGS), op- 
tical spectroscopy both at low and high spectral resolution 
(SMARTS/RC Spectrograph, SMARTS/EcheUe Spectrograph, 
ESO/HARPS), optical polarized spectroscopy (ESPaDOnS), 
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and NIR and optical photometry (SMARTS/ Andicam time se- 
ries). 

Results based on the SMARTS and ESPaDOnS data, includ- 
ing Zeeman-Doppler brightness and magnetic fields maps, are 
presented by Donati et al. ( 2011) ) in a companion paper. In this 
paper we present the analysis of the X-ray ChandrafHETGS ob- 
servation, taking advantage, in the interpretation, of the results 
from the optical data. Magnetic field extrapolation, based on the 
Zeeman-Doppler maps, will be described in forthcoming papers 
and will allow us to derive a three-dimensional description of 
the open (wind-bearing) and closed (X-ray blight) field lines, 
and of the magnetic field structures through which accretion oc- 
curs ( Jardine et al.„200 8), for a more detailed comparison with 
the present X-ray observations. 

We report the main properties of V2129 Oph in Sect. [2] 
Sect. [3] contains the analysis of the X-ray data, including short 
time scales variability (Sect. 3.3 1, and long time scales variability 



in which the analysis of high resolution X-ray spectra is possible 
(Sect. [3!4l l. Our main results are summarized in Sect.[4]and then 
discussed in Sect.|5] 

2. Properties of V2129 Oph 

V2129 Oph is a CTTS belonging to the p Oph star forming re- 
gion at a distance of 120 pc. It is a K5 star, with M - 1.35 M0, 
= 2.1Ro 



R^, = 2.1 Rq, (see D onati et al.|[2()TT 
lower mass companion, separated by 0.6' 
tio ~ 0.09 with respect to the primary ([McCabe et al 



for details). It has a 
'with a /T-band flux ra- 



2006 1: 



using the evolutionary models of |Siess et al.| ( |2000[ l the mass 
of the secondary is ~ 0.1 Mq. Being located at a sky-projected 
distance of ~ 70 AU from V2129 Oph, we can reasonably as- 
sume that it does not significantly affect its magnetosphere and 
accretion process. V212 9 Oph has strong and variable H g emis- 
sion (EW ~ 10 - 14 A, [Bouvier & Appenzeller|l992t . We se- 
lected V2129 Oph for this campaign because it is bright both 
in the optical and the X-ray band, suffers moderate interstellar 
absorption (Ay ~ 0.6, [Donati et al.|2007]l, and because, thanks 



to its rotational period (f - 6.53 d, Shevchenko & Herbst 



1998 Grankin et al. 2008 1 and inclination angle (; ^ 60 deg, 
Donati et al.||201l] irit can be studied with techniques that ex- 
ploit the rotational modulation of surface features. V2129 Oph 
is the first CTTS for which Zeeman-Doppler maps of the surface 
magnetic fields were obtained (Don ati et al.| |2007) from time 
series of spectropolarimetric observations performed in 2005. 
Following this study Jardine et al.| ([2008 ) extrapolated the sur- 
face magnetic field maps and inferred the characteristics of the 
accretion funnels and of the coronal plasma of V2129 Oph at the 
time of the 2005 campaign. No coordinated X-ray observations 
were performed at that time, and indeed, previous to our obser- 
vation, V2129 Oph had been only observed in X-rays by ROSAT 
with low sensitivity and spectral resolution. 

During our 2009 observing campaign we acquired a new 
ESPaDOnS time series during 1-14 July, almost simultaneously 
with the Chandra observation (27-29 June 2009, see Sect. [3|. 
The derived properties of V2129 Oph at the time of the Chandra 
observation are presented by Donati et al. (201 1 1. The presence 
of a hot spot covering ~2.5% of the stellar surface and located 
at high latitude (~ 60 deg) was indicated by Zeeman-Doppler 
maps of the Can IRT excess emission. This excess emission is 
powered by the accretion process, and therefore indicates the lo- 
cation of the footpoint of the main accretion stream in the visible 
part of the stellar surface. Donati et al.ffind that the radial compo- 
nent of the magnetic field has its peak intensity (~ 4kG) in this 
same high latitude region, thus supporting its identification as 



the base of the accretion funnel, i.e. the magnetic flux tube that 
channels the accreting material. This accretion funnel is likely 
trailing the hot spot (i.e. its base), during the stellar rotation, as 
indicated by the periodic red-shifted absorption observed in the 
profiles of the Balmer lines. The intensity of the magnetic field 
suggests that the inner disk is truncated at a distance of 7.27?*. 
The accretion rate of V2129 Oph appears to be quite stable dur- 
ing the observing campaign (apparent variations within a factor 
2) with an average value of log M - -9.2, with M in units of 
Moyr-i. 

The upper part of Fig. [T] shows surface maps of the excess 
emission and of the radial component of the magnetic field from 
[Donati et al. ( 201 1) . The maps are shown as viewed from the 
observer at rotational phases corresponding to the beginning and 
to the end of each of the two Chandra exposures comprising our 
X-ray observation. 



3. Chandra observation and data analysis 

V2129 Oph was observed by Chandra/HETGS on 27-29 June 
2009, for a total exposure of 200.4 ks. The observation, per- 
formed during two consecutive orbits, was divided into two seg- 
ments (obs. id. 9943 and 9944), with durations of 98.9 and 
101.7 ks, and separated by 88.6 ks. The rotational phases of 
the surface maps of excess Can emission and radial magnetic 
field ( [Donati et al.||20lT] l, corresponding to starts and stops of 
the Chandra exposures, are shown in the upper part of Fig. [T 
Rotational phases referenced throughout this paper are referrec 
to the ephemeris defined in [Donati et al.[ ( [2007| l. 



3.1. Data preparation 



The HETGS on board of Chandra ( [Canizares et al.[2005| l is com- 
posed by two gratings, the High Energy and the Medium Energy 
Gratings (HEG and MEG). The HEG and MEG P' order spectra 
cover the 1-15 and 2 - 30 A ranges, with resolutions of 0.012 
and 0.023 A(FWHM), respectively. 

We reprocessed the Chandra/HETGS data using CIAO 4.2 
(CALDB version 4.3.1) starting from the level 1 event file and 
following the standard processing threads described in the ob- 
servatory website. For the analysis we used both the dispersed 
HEG and MEG photons, and the 0* order image. For the 0* or- 
der image, source photons were extracted from a circular region 
centered on the source position (with radius of 5px or ~ 2.5 ", 
expected to encircle > 95% of the source photons), while the 
background was taken from a concentric annulus (with 10 and 
15 px as inner and outer radii). HEG and MEG dispersed spec- 
tra were extracted from the standard rectangular regions aligned 
with the dispersion directions. 

3.2. Zero order image and CCD-resolution spectrum 

We inspected the shape of the point spread function (PSF) of the 
0* order image in order to assess the possible contribution to the 
detected X-rays from the low-mass companion to V2129 Oph. 
With a 0.6" separation, and considering the Chandra spatial res- 
olution, the companion might be marginally resolved if X-ray 
bright. The PSF, however, did not show evidence of asymmetry. 
We will therefore assume that the observed X-rays are all from 
V2129 Oph. This is in agreement with the expected low X-ray 
luminosity of th e secondary that, given the mass-Z^x relation of 
IPreibisch et al.| |2005| l, is ~ 8 x lO^'^ergs^', about 3% of the 
observed Lx- 
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Fig. 1. Chandra light curve in the 0.5-10.0 keV band combining events from the zero-order image and the HEG and MEG first order 
spectra, with 2 ks time bin. Shown above the light curve are images of the excess emission due to accretion (upper row) and of the 
radial component of the magnetic field (lower row) as reconstructed with quasi-simultaneous Zeeman-Doppler imaging by Donati 
et al. ( 201 Colors scales are as in Donati et al. (red and blue in the lower row indicate positive and negative B-fields, respectively). 
The images refer to rotational phases at the beginning and at the end of each of the two X-ray exposures. 



5 

„ 4 

in 
O 

2 

^ 1 






EMg (9.6 MK) 
EMh ( 34 MK) 



50 



100 150 200 

Time [Ks] since JD 2455009.72682 



250 



300 



Fig. 2. Emission measure vs. time. Red and blue points refer to the EM of the cool and hot components, respectively. See Sect, 
for details. 



3.3.2 



We then analyzed the CCD ACIS spectrum obtained from 
the undispersed 0* order image with the main goal to constrain 
the amount of absorption affecting X-rays from V2129 Oph. 
Absorption, parametrized as usual by the equivalent hydrogen 
column density A^h. alters the shape of the X-ray spectrum 
mostly at energy below 1 keV (for modest values of A^h) and, 
in this energy range, the 0* order spectrum has a higher signal 
to noise ratio than HEG and MEG spectra. 

We considered the 0* order spectrum integrated over the 
whole Chandra observing time (we found no evidence of A^h 
variations within the Chandra observation, see Sect. |3.3| l. The 
spectrum, binned between 0.5 and 10 keV so to have a signal 
to noise ratio of 10.0 per bin, was analyzed using XSPEC. A 
two-temperature thermal model (apec, in XSPEC), subject to in- 
terstellar absorption (wabs) was found to adequately fit the ob- 
served spectrum. Elemental abundances were kept frozen to the 



reference values obtained by Maggio et al. ( 2007| l from the X- 
ray emission of young stars of the Orion Nebula Cluster (ONC). 
We evaluated the statistical quality of the fits from the null prob- 
ability implied by the and the number of degrees of freedom. 
The parameters of the emission model are: Ti = 9.6;^[ q, T2 = 
34+^ MK, log£Mi = 52.94 ± 0.11, and logEMz = 53.21+°°^ 
(with EM in units of cm"^). The best fit value for the hydrogen 
column density, A^h, is (8^4) x lO^'^cm"^, fully compatible with 
the value estimated from Ay using a standard conversion rela- 
tion: A^H = 1-6 Av X 10^' ( |Vuong et al.||2003| l. Finally we note 
that the best fit A^h value does not change significantly if the 0* 
order spectrum is fit adopting the abundances that result from 



our analysis of the dispersed spectrum (Sect. 3.4.2 1. 
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3.3. Short term variability 

In order to investigate time-variability on short timescales we 
combined the 0* order undispersed photons with the P' order 
HEG and MEG spectra, collected during selected time-intervals. 
Combining these three datasets allows us to increase the signal 
to noise ratio and hence to explore time scales down to a few 
hours. 



3.3.1. X-ray light curve 

We first constructed a combined background-subtracted light 
curve from the arrival times of photons in source and back- 
ground regions of both the 0* order image and the HEG and 
MEG dispersed spectra. For the 0* order image, source and 
background photons were extracted as indicated in Sect. |3.1 
For the HEG and MEG spectra, in order to increase the signa 
to noise ratio, we selected non standard regions: order source 
events were extracted within 2" -wide strips along the direction 
of dispersion; background events from strips, extending from 7.2 
to 14.4"perpendicular to dispersion direction, on each side of 
each spectrum. Considering the two observing segments, we col- 
lected, after background subtraction, a total of 5852, 2034.5, and 
5341.25 photons (in the E=0.5-10keV band) from the 0* order 
image, the HEG spectrum, and the MEG spectrum, respectively. 

Figure [l] shows the 0* H- P' order light curve of our two ob- 
serving segments, using 2ks binning. The time elapsed since 
the beginning of the observation and the rotational phase of 
V2129 Oph are given in the lower and upper x-axis, respectively. 

The light curve shows a clear flare toward the middle of 
the first segment, lasting ~10 hours and with peak count-rate 
~ 3 times higher than the pre-flare level, and what appears as 
a smaller flare toward the end of the same segment. During the 
second segment a seemingly more gradual variability is observed 
with a slow decay in the first ~35 ksec, and a similarly slow rise 
in the last ~35 ksec. 



3.3.2. Time resolved spectroscopy 

A better insight on the physical origin of the observed variabil- 
ity might come from a time-resolved spectral characterization of 
the X-ray emission. To this end we have extracted spectra in time 
intervals containing a total of 900 counts in the combined 0* or- 
der and P' order (HEG -n MEG) light curve. We oversampled the 
light curve so that each /* interval overlaps with the (/+1)* ones 
but is fully independent from the (/+2)* intervals. For each in- 
terval three spectra were extracted: the 0* order ACIS spectrum, 
and the two l '^' order HEG and MEG spectra. These spectra were 
binned between 0.5 and 10 keV so to have a signal to noise ratio 
of 3.0 per bin. Simultaneous spectral fits to the three spectrc[^in 
the 0.5 - lOkeV band were then performed with XSPEC, using 
both the ;t^^-statistic and the Cash-statistic. 

We adopted 1- and 2- temperatures thermal models (apec, in 
XSPEC), including the effect of interstellar absorption (wabs), 
with plasma abundances fixed to the values of Maggio et aT] 



(2007jl. Fits were performed both leaving the absorption, i.e. the 
hydrogen column density A^h, as a free parameter and fixing A^h 



' Note that for the analysis of the time-average spectrum we used 
the 0* order data only. This was justified by the higher signal to noise 
ratio of the 0* order spectrum at low energies and the need to avoid 
cross-calibration issues. For the lower signal spectra discussed here the 
statistical uncertainties are always larger than the systematic ones and 
cross-calibration issues are therefore less important. 



to 8 X 10^" cm"^, the value obtained in Sect. 3^ for the time av- 
eraged 0* order spectrum. The statistical quality of the fits was 
evaluated by considering the null probability associated to the 
and the number of degrees of freedom. Two thermal com- 
ponents were generally required to obtain acceptable fits, espe- 
cially when fixing the A^h- No evidence of variability of the ab- 
sorption was found and we thus favored models with fixed A^h- 
Although two temperatures were always required for an accept- 
able fit, no evidence for a significant evolution of the tempera- 
tures was found (the only exception being the time interval corre- 
sponding to the rise phase of the flare where a different modeling 
approach indicates a higher temperature for the hot component, 
see Sect. 3.3.3 i. Statistically acceptable fits could be obtained 
by fixing the temperatures to those obtained from the average 
spectrum (Sect. 3.2 1. This simply indicates that the Emission 



Measure Distribution (EMD) of the plasma, while varying sig- 
nificantly in time, can at all times be schematically represented 
by two dominant thermal components at fixed temperatures and 
that no significant additional component can be discerned with 
the limited resolution and signal of the spectra in individual time 
interval^ Fixing the A^h and the two temperatures, the only free 
parameters left in the fits are the emission measures, EM, of 
the two components. Figure |2] shows the time evolution of the 
cool and hot EM. A hardening of the spectrum at times of high 
count-rates is observed (and is also evidenced by a simpler hard- 
ness ratio analysis) both during the flares (see Sect. 3.3.3 1 in the 



first observing segment and in the more gradual evolution in the 
second segment. 

From the above analysis, the absorption-corrected X-ray lu- 
minosity of V2129 Oph ranges, during our Chandra observation, 
from 2.1 to 5.5 X 10^° erg s ' in the 0.5 - lO.OkeV band. 

3.3.3. Flare analysis 

Flares offer the opportunity to constrain the spatial distribution 
of coronal plasma. In particular, the length of the magnetic loops 
where flares occur can be estimated from the evolution of ob- 
servable plasma properties during both the rise and the decay 
phases (Reale et al. 1997; Reale 2007). We have thus analyzed 
in detail the large flare detected during the first observing seg- 
ment. 

In order to use the diagnostic based on the decay phase, the 
heating occurring during the decay must be characterized. This 
is done by studying the time-evolution of the plasma tempera- 
ture and density. In our case, however, a time-resolved spectral 
analysis yielded no evidence of a significant temperature evolu- 
tion during the decay phase, preventing the application of this 
method. 

An estimate of the loop length could instead be obtained 
from the rise phase of the flare (Reale 20071. The quantities used 
to derive the loop length are: the time fo at which the loop tem- 
perature reaches its maximum, the corresponding plasma tem- 
perature To, the time at which the plasma density is maximum, 
fM, and the corresponding plasma temperature, Tm. Both to and 
fM are relative to the beginning of the flare and the temperatures 
refer to the apex of the loop. Two independent estimates of the 
loop half length L are then given by (eq. 12 and 13 of Reale 
120071 ): 



L^3»i'2r,y^M 



2.5 



in¥ 



Tn^ (Im - to) 



(1) 



- We have also repeated the analysis with a coarser time-binning, i.e. 
adopting intervals that contain 2000 counts, instead of 900, with analo- 
gous results. 
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Fig. 3. Upper panel: Chandra X-ray light curve (the same shown 
in Fig.[T]i of the flare. Middle panel: Hardness ratio (HR), defined 
as the ratio between events in the 3.0 - 8.0 keV band and events 
in the 0.5 - 3.0 keV, vs. time. Lower panel: Temperature of the 
hot plasma component during flare evolution. Red points mark 
the maximum temperature and the temperature measured at the 
time in which the count rate is maximum. 



where *P is Tq/Tm, L is in units of lO** cm. To is in units of 10^ K, 
f() and fM are in ks. 

We applied these diagnostics to the flare of V2129 Oph by 
analyzing the 0* -n P' order HEG and MEG events. The times 
to and fM were estimated by inspecting the time-evolution of the 
count-rate and hardness ratio (HR, the ratio between the count 
rate in the hard and in the soft band). The time of maximum den- 
sity corresponds to the time of maximum count rate: this latter 
is indeed approximately proportional to the EM of the plasma in 
the flaring loop and therefore, assuming that the loop volume is 
constant, a proxy of the plasma density. The time of maximum 
temperature corresponds to the time of maximum HR, since the 
HR increases for increasing temperature of the flaring plasma. 

The upper panel of Fig. |3] shows the observed count rate in 
the 0.5 - 8.0 keV X-ray band during the flare. We estimate that 
the flare starts 26 ± 1 ks after the beginning of the Chandra ob- 
servation, and that the density peak occurs at Tm = 43 ± 1 ks. The 
middle panel of Fig.[3]shows the evolution of the HR, evaluated 
adopting 0.5-3.0 and 3.0-8.0keV as the soft and hard bands so 
to maximize the sensitivity to high temperatures. The maximum 
temperature occurs at fo = 35 + 1 ks. The beginning of the flare 
and the times of maximum temperature and density are indicated 
in Fig. [3] with vertical gray bands. 

We determined the temperatures Tq and Tm through time re- 
solved spectral fittings. The adopted time intervals and the re- 
sulting extracted spectra are the same as defined in Sect. |3.3.2| 



Being interested in the temperatures of the flaring plasma we 
must subtract the X-ray emission from the "quiescent" corona. 
We do so by assuming that this latter remains constant during the 
flare and is well represented by the X-ray spectrum registered 
just before the beginning of the flare (in particular, with refer- 
ence to Fig. |2] during the second interval centered at f ~18 ks). 
We thus first derived a spectral model for the quiescent emis- 
sion, and then fitted the spectra of the "flaring" time-intervals as- 
suming a spectral model composed by the sum of the quiescent 
spectrum (with fixed parameters) and an optically-thin isother- 
mal plasma component (apec, in XSPEC) subject to the same 
interstellar extinction as derived from the time-averages spectral 
analysis (A^h = 8 x 10^° cm"^). This latter component, describ- 
ing the flaring plasma, is characterized by two fit-parameters, 
the temperature T and the emission measure EM. The resulting 
best-fit temperatures are shown in the lower panel of Fig. [3] with 
thicker red symbols marking the values of the maximum temper- 
ature, 80^3Q MK, and of the temperature at the time of maximum 
density, 37;^^ MK. Note that these values are averages over the 
real temperature distribution of the flaring plasma, and do not 
directly correspond to the Tq and Tyi used in eq. [T]that are the 
temperatures at the loop apex. Tq and Tm can however be ob- 
tained from the meas ured tempera tures with the transformation 



Reale 



\2mi\, yielding T,, = 2W_^^mi 



formulas reported in 
andTM = 81+24MK. 

From the above estimates of fo, Jm, Tq, and Tm, we finally 
derive two estimates for the length of the flaring loop L with the 
formulas in eq. [T] In order to estimate uncertainties on L, we 
performed Monte Carlo simulation^ of the distribution of the 
input parameters to, tu. To, and Tm, and inferred the L quantile 
corresponding to +lcr. The uncertainties on L are dominated by 
those on To and Tm- In particular the poorly constrained upper 
limit on To yields a poorly constrained upper limit also on L. 
Both expression for L provide very similar results, supporting 
the validity of the analysis. The best values for L are 10-20 Rq, 
with a lower limit at 6 Rq, and a virtually unconstrained upper 
limit. Summarizing we found that the flare occurs in a large-scale 
magnetic loop with half length L~ 10 - 20 R© = 5 - 10/?*. 

We note that there is no evidence of flare occultation by the 
disk or by the star itself because of stellar rotation, and that there 
is no evidence of higher absorption during the flare. This indi- 
cates that during its evolution the flaring loop always remains 
observable and that no significant portion of infalling material 
passes along the line of sight. 



3.4. Dispersed spectra analysis and long term variability 

The dispersed HEG and MEG spectra give access to individual 
line fluxes and allow us to derive detailed information on the 
X-ray emitting plasma such as the emission measure distribu- 
tion iEMD), the element abundances, and the plasma density. 
The EMD is a powerful tool to monitor the amount of plasma 
in small temperature bins and can only be derived from individ- 
ual emission lines, not measurable in low resolution X-ray spec- 
tra. We measured line fluxes and performed our complete spec- 
tral analysis separately for the spectra gathered in each of the 
two observing segments. We can therefore explore, with the de- 
tails offered by high resolution spectroscopy, variability on time 
scales of ~ lOOks. Shorter time scales cannot be investigated 
with these diagnostic tools because of the limited signal to noise 



^ This procedure was preferred to simple error propagation which 
requires uncertainties to be small, since this is not the case for our two 
temperatures. 
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Fig. 4. V2129 Oph X-ray spectra observed with Chandra (HEG+MEG P' orders) during the two observing segments. For clar- 



ity reasons both spectra are rebinned by a factor of 3 (bin size 
wavelengths by 0.05 A. 



0.015 A), and the spectrum of seg. 2 is shifted toward longer 



ratio. We also analyzed the spectrum collected during the entire 
observation so to obtain time-averaged plasma properties with 
the highest possible accuracy. 

In order to increase the signal to noise ratio we added the 
MEG and HEG I'*' order spectra, having first rebinned the HEG 
spectra to the MEG wavelength grid. With this procedure we 
obtained two dispersed spectra, one for each of the two observ- 
ing segments, totaling 3750 and 4274 net counts, respectively. 
Figure. |4] shows the resulting spectra with labels identifying the 
strongest emission lines. The subsequent analysis was performed 
using the IDL package PINTofALE v2.0 ( |Kashyap & Drake 
2000j ), and the CHIANTI v4.2 atomic database ( |Young et al. 
2003 I. 



3.4.1. Flux measurements 

We measured the fluxes of the strongest emission lines in the ob- 
served X-ray spectra of V2129 Oph. Along with the line fluxes 
we also measured continuum fluxes in selected wavelength in- 
tervals free of significant line contamination (see below). 

Line fluxes were obtained through direct integration of the 
binned spectra. For the integration we considered wavelength in- 
tervals centered around the predicted wavelength of each lin^ 
For each line we subtracted the continuum contribution, esti- 
mated from nearby wavelength intervals free of evident emis- 
sion lines, and corrected the resulting flux for the fraction of the 
line-profile, assumed Gaussian, outside of the integration do- 
main. Effective areas and exposure times were finally used to 
convert the observed line counts into fluxes. The observed net 
fluxes are listed, for each observing segment and for the entire 
observation, in Table [T] For each line we indicate the ion gen- 
erating it; line blends are indicated by listing the different ions 
contributing to the emission at that wavelength. We also list the 
temperature of maximum formation, T^ax, of each measured line 
(for line blends T^^j^ is obtained from the sum of the emissivi- 



* Given the discrete nature of the energy binning and the need to sum 
over an integer number of bins, the full width of the integration inter- 
vals is line-dependent, and ranges between 3 and 5 wavelength bins, 
corresponding to 2.2 and 3.4cr of the line profiles. 



♦seg. 1 
- ♦ seg. 2 



^, 10- 



10" 



Ifi 1 1 



10 

wavelength (A) 



15 



20 



Fig. 5. Comparison of measured continuum and lines fluxes in 
the two observing segments. Only features detected at the 2cr 
level in at least one segment are plotted. 

ties of the blended lines). T^ax ranges between 2 and 16 MK. 
This range of maximum formation temperatures indicates that 
there is a significant amount of plasma at a few MK, that did not 
emerge from the analysi s of the time averaged low resolution 0* 
order spectrum in Sect. 



3.2 



of from the analysis of the and 
heavily rebinned P' order spectra in Sect. |3.3| In addition to the 
strong lines, we also included in Table[T]other well known lines, 
even if their fluxes are compatible with zero; upper limits are 
anyway useful constraints for the EMD reconstruction. 

In addition to line fluxes we also measured continuum fluxes 
in selected wavelength ranges. These measures are useful to con- 
strain: a) the source metallicit}]^ and b) and the emission mea- 
sure of the hottest plasma componentaj We measure the con- 
tinuum emission of V2129 Oph in wavelength intervals that are 



' Line fluxes are proportional to abundances, while the continuum is 
not. The relative strength of line and continuum emission is therefore a 
measure of the plasma abundances. 

^ The major contribution to the continuum emission is due to 
bremsstrahlung, whose emissivity increases for increasing tempera- 
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Table 1. Measured fluxes of emission lines and continuum intervals in the X-ray spectrum of V2129 Oph. 









seg. 1 


seg. 2 


seg. 


1-^2 


A" 


Ion 


log r* 

max 


Flux'" 


Flux'" 


Flux" 


Lines 


5.04 


S XV 


7.20 


1.2 ± 0.7 


0.2 ± 0.6 


0.7 


± 0.4 


6.18 


Si XIV Si XIV 


7.20 


0.8 ± 0.3 


1.3 ± 0.4 


1.1 


± 0.3 


6.65 


Si XIII 


7.00 


0.9 ± 0.4 


1.0 ± 0.4 


1.0 


± 0.3 


6.69 


Si XIII Si XIII 


6.95 


0.0 ± 0.3 


0.1 ± 0.3 


0.0 


± 0.2 


6.74 


Si XIII 


7.00 


0.9 ± 0.4 


1.0 ± 0.3 


0.9 


± 0.2 


7.11 


Mg XII Mg XII 


7.00 


0.3 ± 0.3 


0.7 ± 0.3 


0.5 


± 0.2 


8.42 


Mg XII Mg XII 


7.00 


1.1 ± 0.4 


1.5 ± 0.4 


1.3 


± 0.3 


9.17 


Mg XI 


6.80 


0.9 ± 0.5 


0.7 ± 0.5 


0.8 


± 0.3 


9.23 


Mgxi 


6.80 


0.5 ± 0.5 


0.0 ± 0.3 


0.1 


± 0.3 


9.31 


Mg XI 


6.80 


0.0 ± 0.4 


0.3 ± 0.4 


0.2 


± 0.3 


9.48 


Ne X Ne x Fe xxi 


7.00 


0.2 ± 0.4 


0.8 ± 0.5 


0.5 


± 0.3 


9.71 


Ne x Ne x 


6.80 


1.2 ± 0.6 


0.9 ± 0.5 


1.0 


± 0.4 


10.24 


Ne x Ne X 


6.80 


2.5 ± 0.7 


2.7 ± 0.7 


2.6 


± 0.5 


11.74 


Fe XXIII 


7.20 


0.0 ± 0.6 


0.8 ± 0.7 


0.4 


± 0.4 


11.77 


Fe XXII 


7.10 


1.6 ± 0.9 


0.6 ± 0.7 


1.1 


± 0.5 


12.13 


Ne X Ne x 


6.80 


19.9 ± 2.4 


15.7 ± 2.1 


17.8 


± 1.5 


12.26 


Fe XVII 


6.80 


0.5 ± 0.8 


0.9 ± 0.8 


0.7 


± 0.5 


12.28 


Fe XXI 


7.00 


1.1 ± 0.9 


1.6 ± 0.9 


1.4 


± 0.6 


12.83 


Fe XX Fe xx 


7.00 


1.3 ± 1.2 


1.4 ± 1.2 


1.3 


±0.8 


13.45 


Neix 


6.60 


8.1 ± 2.4 


3.6 ± 1.8 


5.8 


± 1.4 


13.55 


Neix 


6.55 


2.3 ± 1.7 


0.3 ± 1.2 


1.3 


± 0.9 


13.70 


Neix 


6.60 


2.8 ± 1.8 


2.8 ± 1.7 


2.8 


± 1.1 


14.20 


Fe xvm Fe xvm 


6.90 


2.2 ± 2.0 


0.6 ± 1.6 


1.4 


± 1.1 


14.37 


Fe xvm Fe xvm Fe xvm 


6.90 


1.2 ± 1.6 


1.9 ± 1.7 


1.6 


± 1.0 


15.02 


Fe XVII 


6.75 


5.5 ± 2.8 


3.7 ± 2.4 


4.6 


± 1.7 


16.01 


Fe xvm O viii viii 


6.50 


3.9 ±3.1 


7.0 ± 3.5 


5.5 


± 2.1 


16.78 


Fe XVII 


6.70 


4.8 ± 3.7 


2.3 ± 3.1 


3.5 


± 2.1 


17.05 


Fe XVII 


6.70 


3.0 ± 3.7 


5.0 ±4.1 


4.0 


± 2.4 


17.10 


Fe XVII 


6.70 


1.7 ± 3.4 


3.8 ± 3.9 


2.8 


± 2.2 


18.63 


Ovii 


6.35 


1.5 ±4.5 


3.4 ± 5.0 


2.4 


± 2.8 


18.97 


VIII O VIII 


6.50 


45.3 ± 13.2 


18.3 ± 9.3 


31.6 


±7.5 


Continuum 


[2.01-^2.51] 






4.9 ± 1.0 


6.6 ± 1.1 


5.8 


± 0.7 


[2.51^2.99] 






5.2 ±0.8 


5.7 ± 0.9 


5.5 


±0.6 


[3.38^3.64] 






3.8 ± 0.7 


5.7 ± 0.8 


4.8 


±0.5 


[3.78^3.88] 






2.7 ± 0.6 


2.8 ± 0.6 


2.7 


± 0.4 


[4.33^4.67] 






5.6 ± 1.1 


8.7 ± 1.3 


7.2 


±0.8 


[6.29^6.53] 






5.4 ± 0.6 


6.5 ± 0.6 


5.9 


± 0.4 



" Wavelengths (A). * Temperature (K) of maximum emissivity. '^ Observed fluxes (10 * ph s ' cm ^) with uncertainties at the 68% confidence level. 



free of significant line contributions (in particular we selected in- 
tervals in which the line contribution with respect to the contin- 
uum expected is below 10%, assuming solar metallicities and for 
any temperature ranging between 10'* to 10** K). The continuum 
emission in these wavelength intervals was estimated by inte- 
grating the observed counts, and by subtracting the background 
contribution. The selected intervals and the corresponding mea- 
sured fluxes are listed in Table [T] 

Figure |5] shows a comparison of the measured fluxes (both 
lines and continuum) in the two observing segments. We only 
plot, as a function of wavelength, the fluxes of features detected 
at the 2cr level in at least one segment. Fluxes at short wave- 
lengths tend to be, on average, higher in seg. 2 than in seg. 1, 
while the opposite results for fluxes at long wavelengths (these 
differences are also marginally observable in the plot of the two 
spectra, Fig.|4]i. In many cases this difference is significant at the 
Icr level. The analysis of the low resolution spectra has excluded 



tures. The continuum emission level is therefore highly sensitive to the 
hottest plasma components. 



significant variability of the absorption, hence we can reasonably 
ascribe differences in line ffuxes to differences in the EMD. 

The two coolest lines detected at the 2cr level in at least one 
segment are the Oviii Lya line at 18.97 A and the Neix reso- 
nance line at 13.45 A, mostly emitted by plasma at 3 - 4MK. 
These lines show flux differences between the two segments of 
1.7 and 1.5cr, with seg. 1 fluxes higher than seg. 2. Combining 
the two differences we obtain a very low probability (0.3%) that 
the EM of the plasma at 3 - 4 MK producing this lines has not 
changed. Continuum fluxes, that probe variations of the hottest 
plasma components, are systematically higher during seg. 2 than 
during seg. 1, firmly indicating that the EM at T > lOMK is 
higher during seg. 2. 

3.4.2. Emission IVleasure Distribution 

We derived the EMD and abundances of the plasma responsible 
for the X-ray emission of V2129 Oph by applying the Markov- 
Chain Monte Carlo method of [Kashyap (fe Drake] ( |1998| l to the 
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Table 2. EMD and abundances of V2129 Oph. 
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Fig. 7. EMD obtained from the flux measurements of the two 
observing segments. 



measured continuum and line fluxes in Table[T] Measured fluxes 
were converted into unabsorbed fluxes adopting the hydrogen 
column density A^h of 8 x 10^" cm"^ derived from the 0* order 
spectrum (see Sect. 3.2 1. For the EMD and abundance recon- 



struction we considered all the measured fluxes discarding only 
the lines that depend on plasma density, i.e. the intercombination 
and forbidden lines of He-like triplets. 

Starting from the flux measurements and considering the 
peak temperature at which each line is produced (see Table[T]i we 
reconstructed the EMD over a regular logarithmic temperature 
grid with 0.15 dex bins in the log T range between 6.2 and 7.5. 
Abundances were estimated for those elements for which at least 
one line was measured. Note that the reconstructed EM values 
associated to bins at the boundary of the chosen temperature do- 
main may be overestimated: artificially limiting the temperature 
range of the EMD reconstruction forces to zero the EMD out- 
side this range. Plasma at higher and lower temperatures may 
however contribute to some of the measured fluxes, and thus 
causing the EM values at the boundary bins to be overestimated 
in order to account for the observed fluxes. 

The EMD and abundances derived from the whole observa- 
tion are reported in Table|2]and the EMD is also plotted in Fig. 6] 
Figure |6] also shows the 2 temperature model obtained from the 
0* order spectrum of the whole observation (see Sect. |3.2|. We 
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"•^■-'-0.04 



Abundances are in solar units ( jAsplund et al.|2005) . EMD values are 
in cm"^ and have been derived assuming that V2129 Oph is located at 
120 pc. 



note that the two components of this 2-T model agree very well 
with the two hottest peaks of the EMD. 

We then derived EMD and abundances of the emitting 
plasma considering separately the two observing segments. For 
both segments unabsorbed fluxes were evaluated assuming A^h = 
because no evidence of absorption variability was 



8xl02°cm-2 

found (see Sect. |3.3.2| ). No significant variation in plasma abun- 
dances can be discerned between seg. 1 and 2, indicating that 
the observed variations of line fluxes can be ascribed to EMD 
variations only. We therefore repeated the EMD reconstruction 
for seg. 1 and 2 freezing the plasma abundances to the values of 
the seg. 1+2 model, and leaving as free parameters the EM val- 
ues. The thus derived EMD for the two segments are reported in 
Table [2j and plotted in Fig. [7] As a consequence of the observed 
systematic differences in measured line fluxes, a difference is 
seen in the cool plasma component] with EMDi > EMD2 at the 
1.5o- level in the bin centered on log T = 6.55 {T ~ 3 - 4MK). 
Note that the significance of the time variability is lower for the 
EMD than for the line fluxes (see Sect. |3.4.1| l because the EMD 
reconstruction procedure involves further uncertainty sources 
(Kashya p & Drake] 1998[ l. However the high confidence level 
inferred for the variation of the flux of cool lines indicates that 
the EMD variability is indeed real. 



3.4.3. Electron density 

The plasma density A^e can be constrained, from the high reso- 
lution X-ray spectroscopy, by analyzing the He-like line triplets, 
and, in particular, the flux ratio between the intercombination 
line, /, and the forbidden line, / ( jGabriel & Jordan] |1969| l. 
Increasing plasma densities result in decreasing //; ratios. The 
intensities of the / and / lines depend on plasma density because 



' Note that the corresponding models obtained without fixed abun- 
dances, have EMD values almost identical to those obtained by fixing 
them, but the associated EMD uncertainties are larger and the differ- 
ences at low temperatures are not significant. 
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Fig. 8. Spectra of the two observing segments in the Ne ix triplet 
region. Observed spectra are rebinned by a factor 2 (bin size - 
0.01 A). Model spectra are obtained using line and continuum 
measurements (see Sect. 3.4.1 1 and assuming a gaussian line 
spread function. 




Fig. 9. Forbidden / to intercombination / line ratio of the Ne ix 
triplet vs electron density. Black curve represents the /// pre- 
dicted ratio. Observed // / values of segments 1 and 2, and the 
corresponding density ranges, are indicated in blue and red re- 
spectively. 



these lines originate from radiative decay of metastable states 
(long lifetime states), whose depopulation is favored by increas- 
ing collision rate (i.e. increasing density), while the resonance 
lines originate from excited states depopulated only by photon 
emission processes. 

Among the He-like triplets in the wavelength range covered 
by the ChandrafHETGS instrument, we identified and measured 
the He-like line triplets of Sixiii (at 6.65, 6.69, 6.74 A), Mgxi 
(at 9.17, 9.23, 9.31 A), andNeix(at 13.45, 13.55, 13.70A).The 
analysis of these triplets can constrain the electron density of 
plasma at the formation temperatures of these ions, ~ 10, ~ 6, 
and ~ 4 MK, respectively. Measured line fluxes for these triplets 
are reported in Table [T] 

The Si XIII triplet has, in all of our spectra, / line fluxes com- 
patible with zero within Icr, providing therefore only a lower 
limit for the /// ratio. We therefore infer that plasma at T ~ 
10 MK has a density lower than ~ lO'^ cm"^, during both seg. 1 
and seg. 2. The Mg xi triplet, mostly emitted by plasma at 6 MK, 
has both / and / line fluxes compatible with zero in all the spec- 
tra, thus providing no constraint on plasma density. 

In Fig. [8] we show the lines of the Ne ix triplet correspond- 
ing to the two observing segments. These lines are produced by 
plasma at ~ 4 MK. In the plot we also show noise-free spectra 
reconstructed on the basis of the measured continuum and line 
counts. The overall flux of the triplet is clearly lower in the sec- 
ond observing segment than in first (see also Table[T]i. 

To obtain an accurate estimation of the /// ratio uncertain- 
ties we considered the measured counts of / and / lines, as well 
as those used to determine the continuum level. Assuming that 
these quantities follow Poisson distributions we derived through 
Monte Carlo simulations the distribution of the /// flux ratio, 
and hence inferred the /// mean value and the quantile cor- 
responding to ±lcr. With this procedure, and having also cor- 
rected for interstellar absorption, we obtained /// = 1 -2;^,' ^ and 
/// > 2.3 for seg. 1 and 2, respectively. The comparison of the 
observed /// values with the predicted ///-ratio vs. plasma den- 
sity functiorj^is shown in Fig. [9] From this relation we derive: 



logA^e = l^.l^'j'g and logA^e < 11-5 for seg. 1 and 2, respec- 
tively, with A^e in cm"-*. 

4. Main results 

The Chandra observation of V2129 Oph has revealed that it 
is an X-ray bright source. The average X-ray luminosity, de- 
rived from our mean EMD, is Lx - 2.7 x lO-'^ergs"' in the 
0.5 - 10.0 keV band. This Lv, combined with the bolometric 



luminosity of V2129 Oph (logLboi/i^o = 0.15, Donati et al. 



2011), yields logLx/^boi = -3.3, a value that puts V2129 Oph 
among the X-ray brightest accreting T Tauri stars when com- 
pared to CTTS in the ONC ( [Preibisch et"aL][2005) . The X-ray 
emission from V2129 Oph is produced by significant amounts of 
plasma spanning a broad temperature range, at least from 2 MK 
up to ~ 30 MK, as evidenced by the reconstructed EMD (see 
Fig. |6]l. In particular the plasma at a few MK, monitored dur- 
ing the first half of the Chandra observation, shows high den- 
sity, logA^e ~ 12.1, consistent with results for almost all the 
other CTTS with measured densities, supporting the idea that 
the presence of high density cool plasma is directly linked to the 
accretion phenomenon. 

4.1. X-ray variability 

The X-ray emission of V2129 Oph is variable on short time 
scales (hours). In the 200 ks Chandra observation we detected 
a clear flare, ~10h long, causing an increase in the X-ray lumi- 
nosity by a factor ~ 3. This flare occurred in a long loop, with 
an analysis of the rise phase giving a best value of the loop half 
length of ~ 5 - IQR*, with a lower limit of 3 R*, and an uncon- 
strained upper limit. 

Beside this clear flare, all the observation is characterized by 
variability similar to that typically produced by coronal activ- 
ity, i.e. with clearly correlated X-ray emission level and spectral 
hardness. This short term variability is in fact mainly produced 
by the hottest (r ^ 34 MK) plasma components. 



^ The predicted ratios, other than on density, also depend, very 
weakly, on the plasma EMD. In deriving the //(-ratio vs. density re- 
lation we adopted the mean EMD for whole observation. Using the 



individual EMD derived for the two segments, however, produces neg- 
ligible differences in the inferred log A'e, smaller than 0.01 dex. 
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The analysis of the high resolution spectra collected during 
the two segments allowed us to search for variability of the X-ray 
emission observed at two different phases during stellar rotation 
and, therefore, emerging from either different magnetospheric 
structures (corona and accretion shock) or from the same struc- 
tures seen with different viewing geometries. We stress that the 
analysis of the high resolution spectra allows us to probe the 
coolest plasma components. This comparison yielded two main 
results. 

1 . Cool plasma EM variability. Lines at long wavelengths were 
significantly enhanced during seg. 1 with respect to seg. 2. 
This finding is apparently not related to variation in the ab- 
sorbing column since no such variability is derived from the 
XSPEC analysis of the heavily rebinned 0* + 1* order spec- 
tra. The EMD of the two observing segments indeed in- 
dicates that this diff'erence originates from a difference in 
the amount of observable emission measure at low temper- 
atures. In particular we observe a difference in the EM at 
r ^ 3 - 4MK (log T = 6.55), with EMDx being higher than 
EMD2, at the 1.5cr level, by more than a factor of 10. The 
difference in the amount of emission measure at a few MK 
is likely not related to flaring activity, since, as indicated by 
the short time-scale analysis, flares contribute to the EMD at 
much hotter temperatures than 4MK (see Fig. [2]). 

2. Plasma density variability. During the two observing seg- 
ments we found different constraints on the density of the 
plasma component at a few MK. In segment 1 the cool 
plasma density can be constrained between 1.3 x 10'' and 
4 X 10'^ cm"-*, while during seg. 2 the density diagnostic is 
compatible with a low density plasma, A^e < 3 x lO" cm""*. 
Although the statistical significance of this difference is low, 
V2129 Oph, is the first CTTS where apparent variations of 
the plasma density over time scales of ~ lOOks have been 
observed. 



5. Discussion 

We now discuss the implications of our observational results on 
the X-ray emission of V2129 Oph, taking advantage of the con- 
strains on the properties and geometry of accretion shock region 
and on the stellar magnetic field, obtained by us simultaneously 
with the X-ray observation ( |Donati et al.|2011 1. 

The analysis of X-ray data allows us to derive the main prop- 
erties of the coronal plasma and of the plasma heated in accre- 
tion shocks, as discussed below. These results will also provide 
constraints for future works dedicated to magnetic field extrap- 
olation, aimed at modeling the coronal plasma and the accretion 
stream structures. 



5.1. Coronal emission 

The broad band X-ray luminosity of V2129 Oph is dominated 
by emission from hot plasma components, most likely of coro- 
nal origin. In 1991 flie ROSAT satelHte detected V2129 Oph 
Casanova et al. 1995) with an X-ray luminosity of 2.5 x 



0^" erg s~\ almost identical to that measured with Chandra in 
2009. V2129 Oph thus appears to have a stable X-ray emission, 
hence coronal activity level, over time scales of ~ 20 yr. 

V2129 Oph shows X-ray variability over short time scales 
(Fig.[T]i. This variability, being mainly due very hot plasmas, can 
be attributed to coronal activity. The X-ray light curve of the first 
observing segment is characterized by an obvious flare. We con- 
strained the half length of the flaring loop, deriving 5 - 10/?* 



as best value (with a lower limit of 3/?* and with an uncon- 
strained upper limit). Considering its length and the inner disk 
radius (~ 7.2 /?*), we cannot exclude the possibility that this loop 
has one footpoint anchored to the inner circumstellar disk. Such 
flaring coronal structures connecting star and inner disk were 
already suggested in young stars of the Orion Nebula Cluster 
( |Favata et al.|2005[ [Getman et al.f 2008), although we note that 
the presence of a disk is not a prerequisite for such large flaring 
loops to exist. They are also found to occur on stars which show 
no evidence for disks dGet man et al.|2008]|Aarnio et al.|20T0l l. 

The second segment of the Chandra observation shows a 
trend in the first part of the segment that might be interpreted 
as the decay phase of a flare, followed in the second part by an 
increase of the count rate and of the hardness ratio that might be 
associated the rise phase of another flare. 

Alternatively the variability might be due to modulation of 
some emission feature on the stellar surface as the star rotates. 
We note that the observed dip in the light curve coincides with 
the passage along the line of sight of the large region of strong 
positive radial magnetic field (in red in Fig[T]l, which is also co- 
incident with an extended dark photospheric spot (not shown, cf. 
Don ati et al.|20 1 1 ). Although we can only speculate that the dip 
in the X-ray light curve is associated with the passage of this 
region, it is not unconceivable that the X-ray emitting corona, 
assumed to be confined by relatively small-scale closed mag- 
netic loops, is disrupted in this region of strong magnetic field 
likely connecting the star with the inner disk (as it is certainly 
true for the field in the accretion spot that lies within the region). 
Figure |2] might further indicate that while this region removes a 
significant fraction of visible hot (34 MK) plasma, it leaves the 
cool (9.6 MK) plasma relatively unaffected thus indicating that 
the two are not co-spatial. 

5.2. Higii density cool plasma and accretion shock(s) 

We found that, during the first segment of our Chandra observa- 
tion, the plasma at a few MK has a high density (log A^e - 12.1). 
Therefore V2129 Oph is a further case confirming that in CTTS 
cool plasma has high density. This plasma component is in- 
trinsically different from the coronal plasma observed in non- 
accreting active stars. The commonly accepted model suggests 
that this plasma component is material heated in the accretion 
shock. In fact, the high infall velocity of the accreting material 
naturally justifies material at such temperatures, of a few MK, 
and densities. 

Plasma heated in accretion shocks is located at the base of 
the accretion stream, just behind the shock surface. Models pre- 
dict that this accreted material, passing through the shock sur- 
face, strongly reduces its infall velocity, but continues its down- 
ward motion across the stellar atmosphere, meanwhile coohng 
down and radiating X-rays ( [Giinther et al. 2007 Sacco et al. 
[2008; Orlando et al. 201^1. 

The association of the cool plasma component at ~ 3 - 4MK 
of V2129 Oph with the accretion shock is greatly reinforced if 
we consider that this temperature is almost exactly what we ac- 
tually expect for the post-shock plasma. The measured strength 
of the magnetic field of V2129 Oph implies that inner disk is 
truncated at LIR^, ( |Donati et aL]|2011| l. A free fall from that 
distance implies a pre-shock velocity of 460 km s"', and, as a 
consequence, a post-shock temperature of 3 MK. 

We will therefore assume that X-rays from V2129 Oph are 
emitted by two plasma components: a) the high density plasma 
component at 3 - 4 MK, heated in the accretion shock and hence 
located at the base of the accretion stream, and b) a low density 
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corona, with T ranging from 2 up to a ~ 30 MK. We will try 
to interpret the time variability of the EM at ~ 3-4 MK and, 
possibly, of the measured plasma density that we observed be- 
tween our two observing segments in terms of different viewing 
geometry of the accretion shock and funnel due to stellar rota- 
tion. Figure[TO]displays a cartoon describing the system and the 
viewing geometry during the two observing segments. 

From a theoretical point of view, it has been debated whether 
X-rays emitted by shock-heated plasma in CTTS are able to 



escape from the shock region, and hence be detected (Drake 
2005 Sacco et"aL]|2010 i. Their escape probability depends on 



the height of the shock surface on the stellar atmosphere and on 
the direction of emission. In fact. X-rays emitted in the direction 
of the accretion stream could be absorbed by the pre-shock in- 
falling material, depending on the accretion stream density and 
dimensions along the line of sight; X-rays emitted perpendicu- 
larly to the accretion funnel could be absorbed by the surround- 
ing stellar atmosphere, depending on the sinking of the hot post 
shock in the stellar atmosphere ( Drake|2005 Sacco et al.|2010[ l. 
The amount of X-rays produced by the accreton shock able to 
escape in the different directions may be estimated through de- 
tailed 3D modeling of the shock region, which is beyond the 
scope of this paper 

In the case of V2129 Oph, we can assume that the post-shock 
plasma is co-spatial with the the optical hot spot, whose loca- 
tion we have determined through Zeeman-Doppler imaging. As 
shown in Fig. [T] and outlined in the cartoon of Fig. 10 this re- 
gion remains visible during both Chandra observing segments, 
but with distinctly different viewing geometries: near the stellar 
limb during seg. 1 and nearer the center of the stellar disk during 
seg. 2. 

During seg. 1 the accretion funnel is observed sideways. In 
this configuration. X-rays produced by the shock heated plasma 
can easily be observed because they can escape from the lat- 
eral part of the funnel. In fact the accretion shock characteristics 
of V2129 Oph, i.e. a pre-shock velocity of 460 kms ', a low 
metallicity, and a pre-shock density of ~ 3 x 10" cm"^ (since 
the accretion shock is expected to be a strong shock, the density 
of the pre-shock material should be one fourth of the measured 
post-shock density), imply that the hot post-shock region is high 
enough on the stellar surface to avoid significant absorption from 
the surrounding chromosphere ( Sacco et al. 2010[ l. Of course X- 
rays from coronal plasma are also observable. The emerging X- 
ray spectrum and the resulting EMD are therefore the superpo- 
sition of the two contributions, with the post-shock component 
dominating the EMD at 3 - 4 MK, and the corona dominating at 
higher temperatures. The observed X-ray spectrum of seg. 1 is 
in fact compatible with an overall low absorption, and with high 
density for the plasma at ~ 4 MK. 

Conversely, during seg. 2 the hot spot almost faces the ob- 
server. Pre-shock infalling material is thus located between us 
and the post-shock region at the base of the accretion funnel 
(see Fig.[TO|. The X-rays emitted in our direction by post-shock 
plasma thus enter the pre-shock accretion stream. Even consid- 
ering that our line of sight intersects a small fraction of the ac- 
cretion funnel, i.e. 15% of stellar radius, the resulting hydrogen 
column density is ~ 7 x 10^' cm"^ (having assumed a pre-shock 
density of 3 x 10" cm""*). Such a large column density would 
absorb 90% of X-rays at ~ 13.5 A (the wavelength of the Ne ix 
triplet). In this viewing geometry essentially no X-ray emitted 
by the shock heated plasma can be observed. Therefore during 
seg. 2 only coronal X-rays can be observed, because they are not 
likely to be blocked by the pre-shock material. The emerging 
X-ray spectrum is therefore almost completely due to coronal 



plasma. The reconstructed EMD of seg. 2, hence misses the high 
EM values at 3 - 4MK, and the the seg. 2 Neix triplet reveals 
low density plasma, compatible with a coronal origin. In this ge- 
ometry the emerging X-ray spectrum suffers again the same low 
absorption of seg. 1. 

We notice that the lowest temperature bin of the EMD (at 
2 MK), even if poorly constrained, appears constant between the 
two segments. This does not contradict the above picture for two 
reasons: 1) the error bars on the EM might hide even a large 
variation, 2) HD simulations predict that accretion shocks onto 
CTTS produce sharply peaked EMD ( |Sacco et al.1[20T0l l that 
might not contribute significantly to this EM bin. 

The proposed scenario is able to explain the observed varia- 
tions of the EMD and density with significant variations of the 
extinction suffered by the X-rays from the accretion shock. Note 
that this is compatible with the lack of variability in the absorp- 
tion measured from the X-ray data. In fact an increase of the 
absorption in the X-rays of the post-shock plasma only cannot 
be detected because of the presence of the underlying weakly 
absorbed X-ray coronal emission. 

Note that in the optical band, the accretion indicators (e.g. 
the Ca n and He i lines) reach their maximum at phases 0.55-0.7, 
corresponding to rotational phases at which the hot spot is di- 
rected toward us. We have here observed that the opposite hap- 
pens for the accretion-driven X-ray emission: it is maximum 
when the hot spot is near the stellar limb and minimum when 
the hot-spot is facing us. We can make sense of these two op- 
posite behaviors because the infalling material should be almost 
completely in the gas phase. Gas is indeed transparent to opti- 
cal radiation, so that the Can and Hei emission can reach us 
through the funnel, but can easily absorb soft X-rays with ener- 
gies £■ < 1 keV, such as those produced in the accretion shock. 

We have argued that the plasma at 3 - 4 MK observed on 
V2129 Oph is likely to be shock-heated material. The high den- 
sity, the time variability of its EM in agreement with this pic- 
ture, and the coincidence of the temperature with the expecta- 
tion for the post-shock material reinforce this interpretation. In 
this likely hypothesis we can infer the properties of the shock 
heated plasma, and in particular the accretion rate, needed to ex- 
plain the observed soft X-ray flux. From the observed flux of the 
Ne IX resonance line measured during seg. 1 , and assuming a pre- 
shock velocity of 460 km s"' (compatible with a free fall from 
the estimated inner disk position), we derive an accretion rate of 
7 X 10"" M0 yr"'. This is lower by a factor 10 with respect to 
the accretion rate d erived from optical lines, 6 x 10"'"MQyr"' 
( Donati et al. 2011J. This underestimation of M derived from 
X-rays is common for CTTS observed with high resolution X- 
ray spectroscopy ( Argiroffi et al.||2009 Curran et al.| |201l]l. As 
suggested by iGiinther et al. (2007i and investigated by (Sacco 
|et al.||2010| , this discrepancy could be explained assuming that 
the accretion stream is not uniform in density and/or in velocity, 
and that therefore the portion of the accretion stream produc- 
ing observable X-rays is a small fraction of the whole stream. 
Moreover Curran et al. ( 2011) , comparing mass accretion rates 
of CTTS derived from optical diagnostics with those derived 
from high resolution X-ray spectra, argued that, even if X-rays 
provide significantly underestimated values, a weak correlation 
between the two quantities may exist. The values that we obtain 
for V2129 Oph support this tentative correlation. 

Finally, having argued that, at least during seg. 1, the EMD 
bins at low temperature are due to the accretion-shock and those 
at high temperature to coronal plasma, we estimate separately 
the X-ray luminosities of the two components. We do so by 
assuming that the EMD below and above log T - 6.8 are en- 
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Fig. 10. Upper panels: Cartoon describing the spatial distribution of the two X-ray emitting plasma components of V2129 Oph. Left 
and right cartoons correspond to the two viewing angles of the two Chandra observing segments, with observer in the rightward 
direction. Red region marks the post-shock high density plasma at 3 - 4MK, blue regions indicate low density coronal plasma with 
r ranging from 2 up to ~ 30 MK. Lower panels: EMD corresponding to the two observing segments (see Sect. 3.4.2 1. Red EM bins 
symbolize observed EM values ascribed to post-shock plasma, while blue bins those accounting for coronal plasma. During seg. 1 
the pre-shock material does not block the view of both post-shock and coronal plasma: the EMD is therefore the superposition of 
the EMD two components. During seg. 2 the pre-shock material almost completely absorbs the X-rays of the post-shock plasma 
emitted toward the observer, while coronal emission is mostly unaffected: in this case all the X-rays detected are those produced 
by coronal plasma; the reconstructed EMD, being only those of coronal plasma, hence misses the high EM values at 3 - 4MK 
(reported for comparison with a dashed red line in the EMD plot). 



tirely due, respectively, to shock heated material and to coronal 
plasma. From the EMD of seg. 1, the resulting Lx are 3.4 x 10^^ 
and 2.2 x lO^'^ergs for the shock heated and the coronal 
plasma, respectively. Even when we do observe X-rays from 
shock heated material (seg. 1), the coronal X-ray luminosity is 
thus ~ 90% of the total Lx. 

A similar conclusion can be reached by comparing the EMD 
reconstructed from individual line and continuum fluxes, with 
the 2 temperature model obtained from the spectral fitting of 
the CCD-resolution ACIS spectrum (see Fig. |6]l. The temper- 
atures in this latter model correspond to the two hottest peaks 
of the EMD distribution, both ascribable to coronal plasma. 
This is a further confirmation that X-ray spectra of CTTS in the 
0.5 - 10 ke V band are dominated by coronal emission. The shock 
heated material is probed only by the softest part of the spectrum, 
or by emission lines produced by material at temperatures of few 



MK, that contribute little flux to the overall low-resolution spec- 
tra. 



5.3. V2129 0ph vs. TW Hya 

TW Hya is considered the prototype of CTTS with high-density 
X-ray emitting plasma. TW Hya is certainly not the most repre- 
sentative star of its class, being among the oldest CTTS. Because 
of its proximity, however, it is one of the best studied, especially 
in the X-ray band ([Kastner et al.|2002[|Stelzer & Schmitt^200"4t 



IBrickhouse et al.|2010|lT 

We show in Fig. Ill] the comparison between the differential 



emission measurd^ (t/em) of V2129 Oph (relative to seg. 1) and 



that of TW Hya derived by |Brickhouse et al.| ( |20T0l l. 



' The dem is defined as EMD/ A log T. The use of the dem instead of 
the EMD eases the comparison between models defined over tempera- 
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Fig. 11. Differential emission measure (dem) of V2129 Oph, 
compared to that of TW Hya derived by Brickhouse et al. (2010 1 



The X-ray emission of TW Hya is dominated by plasma at 
2.5 MK, as evidenced by the sharp peak in the dem. This plasma 
component has very high density and is thus material heated in 
the accretion shock. The high temperature tail (T > 10 MK) of 
the dem of TW Hya indicates instead that the amount of coronal 
plasma is significantly lower. 

The thermal structure of the plasma on V2129 Oph dif- 
fers significantly from that of TW Hya. The cool plasma 
on V2129 Oph has much lower EM than on TW Hya and 
V2129 Oph thus hosts a significantly smaller amount of high 
density cool plasma related to the accretion process. Conversely 
the corona of V2129 Oph is significantly more active than that of 
TW Hya, bearing a larger amount of emission measure at high 
temperatures. Other than differing in absolute values, we note 
that V2129 Oph and TW Hya have very different relative contri- 
butions of accretion-driven vs. coronal plasmas. 

We note finally that the pole-on inclination under which 
TW Hya is viewed (/ x 7 deg, |Qi et al.|2004| l raises some ques- 
tions on the observability of accretion-driven X-rays in terms of 
viewing geometry of the accretion shock. In fact, postulating that 
also TW Hya has a hot spot latitude of ~ 60 deg, then its accre- 
tion shock should be viewed at the same angle as the accretion 
shock of V2129 Oph during seg. 2. However accretion-driven X- 
rays from TW Hya are always observed. This difference could be 
reconciled assuming for TW Hya a lower latitude of the accre- 
tion shock (as suggested by Batalha et al.|2002) ). Further obser- 
vational constraints on the magnetospheric accretion geometry 
of TW Hya are however needed to address this issue. 

6. Summary and conclusions 

Our 200 ks Chandra observation proved that V2129 Oph is one 
of the X-ray brightest accreting T Tauri stars, with Lx - 2.1 x 
lO^^ergs^' and logLx/^^boi = -3.3. The hottest components of 
the X-ray emitting plasma showed flaring activity typical of ac- 
tive coronae. Conversely the plasma at temperature of 3 - 4 MK 
evidenced a high density, similar to that observed in almost all 
CTTS with measured density, and contrasting with a coronal ori- 
gin. Moreover the Chandra observation of V2129 Oph for the 
first time gave a hint of time variable density: during the first 
half of the observation the plasma at 3 - 4 MK had high density 



ture grids with different bin sizes since the EM in a given temperature 
range can be simply evaluated as the area under the dem curve. 



(logA^e ~ 12.1), while during the second half it showed lower 
density, although the statistical significance of this density dif- 
ference is marginal. The fluxes of strongest emission lines pro- 
duced by this cool plasma component, that probe the EM, are 
correlated with plasma density: in the first half of the observa- 
tion they were significantly higher than in the second. 

These results confirm that the accretion phenomenon is as- 
sociated with the presence of high density plasma at a few MK, 
which is thus likely to be material heated in the accretion shock. 
They also suggest that variations of the observed density and EM 
of this plasma component can occur over time scales of ~ 100 ks. 

Assuming that this high density plasma is material heated in 
the accretion shock, these observed variations can be interpreted 
in terms of rotational modulation: in fact, the viewing geometry 
of the accretion funnel, as determined through simultaneous op- 
tical observations, changes significantly during the Chandra ob- 
servation. During the first half of the observing time the X-rays 
produced in the accretion shock region are observable, while 
during the second half of the observation the dense portions of 
the accretion funnel near the shock blocks the accretion-driven 
X-rays emitted in the direction of the observer. 
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